Pion transition form factor and distribution amplitudes in large-iV c Regge models 
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We analyze the ir-y*^* amplitude in the framework of radial Regge models in the large- N c limit. 
With the assumption of similarity of the asymptotic Regge p and uj meson spectra we find that 
the pion distribution amplitude is constant in the large- N c limit at the scale Qo where the QCD 
radiative corrections are absent - a result found earlier in a class of chiral quark models. We discuss 
the constraints on the couplings from the anomaly and from the limit of large photon virtualities, 
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I. INTRODUCTION 



The study of exclusive processes in QCD pjj Q has 
been a permanent challenge both on the theoretical as 
well as the experimental side. The production of neu- 
tral pions by two virtual photons provides the simplest 
process where both perturbative and non-perturbative 
aspects of strong interactions can be tested. Indeed, the 
normalization of the corresponding from factor for real 
photons is dictated by the chiral anomaly. In the op- 
posite limit of large photon virtualities the amplitude 
factorizes into power corrections and a soft and scale de- 
pendent distribution amplitude (DA) which gives direct 
information on the quark content of the pion at a given 
scale. Actually, for large momenta the behavior of the 
DA is controlled by the perturbative QCD renormaliza- 
tion group evolution 0, 0, 0] m terms of Gegenbauer 
polynomials as implied by the conformal symmetry (for 
a review see, e.g., Ref. and the high-Q 2 power expan- 
sion. Then at Q 2 — » oo one has the fixed point result, 
e.g., (p( 2 \x,Q 2 ) —* 6#(1 — x). For low scales genuinely 
non-perturbative evolution can be tackled by transverse 
lattice methods [10 (for a review see e.g. Ref. 0). 
The first and second moments of the DA have been com- 
puted on Euclidean lattices fToL Hl|. The QCD sum 
rules have been applied to the leading twist- 2 DA [T^| . 
Measurements of the transition form factor have been 
undertaken by the CELLO 13] and CLEO [3 collabo- 
rations. An analysis of the lowest Gegenbauer moments 



'Electronic address: earriola@ugr.es 

^Electronic address: Wojciech.Broniowski@ifj.edu.pl 



a 2 and a 4 has been carried out by and 
Higher twists have been analyzed in the context of light 
cone sum rules |2fJ . A direct measurement of the DA has 
been presented by the E791 collaboration [2l|. For a con- 
cise review on all these developments see, e.g., Ref. [2^ 
and references therein. 

Most calculations dealing with the pion transition form 
factor and more specifically with the DA involve quarks 
as explicit degrees of freedom. This appears rather natu- 
ral but the principle of the quark-hadron duality suggests 
that it should also be possible to make these calcula- 
tions entirely in terms of the complete set of hadronic 
states. In fact, the large iV c -limit of QCD [23II24I ] makes 
quark-hadron duality manifest at the expense of intro- 
ducing an infinite number of weakly interacting stable 
mesons and glueballs. The large- iV c -limit may in fact 
be regarded as a model-independent formulation of the 
quark model. Actually, chiral quark models are partic- 
ular realizations implementing in a rather natural way 
this large- N c behavior at the one-quark-loop approxima- 
tion and several calculations have been made along these 

lines HE HE I^LEM, 123, [S3, H3, H3, H3 (for a review 

see, e.g., Ref. |3qb Despite the subtleties regar ding the 
correct implementation of chiral Ward identities |35L l8r1| , 
chiral quark models by themselves cannot be "better" 
than the large- N c limit, as they form a particular model 
realization of this limit. Surprisingly, up to now there has 
been remarkable little information on the quark content 
of hadrons based solely and directly on the large- N c limit 
and the quark-hadron duality ideas, without resorting to 
specific low-energy quark models. 

The purpose of this paper is to fill this gap and analyze 
the pion transition form factor in the original spirit of the 
large- N c limit. We impose chiral constraints at low ener- 
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gies and QCD short distance constraints at high energies 
and extract from there the DA in a Regge model with in- 
finitely many resonances where the radial squared mass 
spectrum is assumed to be linear. This near-linearity is 
supported by the experimental analysis of Ref. [33] ■ We 
do the analysis in the absence of radiative corrections 
which must be introduced via QCD evolution, hence our 
calculation corresponds to some reference scale Qq where 
the radiative QCD corrections are absent. Remarkably, 
we obtain the same main result as found previously by us 
in chiral quark models 33, 34], namely, the leading- twist 
pion DA is constant at the scale Qq. Our calculation pro- 
vides an explicit example of quark-hadron duality in an 
exclusive process. 

Up till now the calculations within the framework of 
large--/V c Regge models have been mainly restricted to 
two-point functions [H |H EH El El E! On the 
other hand, large-A c -motivated calculations of three- 
point functions with short distance constraints have been 
carried out with a finite number of resonances, which en- 
abled getting model-independent results for vector meson 
decays [441 |45|, 1461 1471 148| . In this regard there are large 
N c studies of both the pion E^ and the proton [5(j elec- 
tromagnetic form factors. It should also be mentioned 
that, echoing some older ideas of Radyushkin [H]]], the 
light-cone wave functions have recently been computed 
within the holographic approach to QCD based on the 
AdS / CFT correspondence p52l l53| or meromorphization 
ideas [54j. In these works the quark-hadron duality is 
also exploited. 

The outline of the paper is as follows. In Sect. ITT1 we 
review DAs in the large- N c context and fix our notation 
for the remainder of the paper. In Scct. lllll we undertake 
the analysis of the large--/V c Regge model and its con- 
sequences for the pion transition form factor and pion 
DA. Section HVl deals with aspects of the QCD evolution, 
which is a crucial ingredient of analyses of this kind. For 
completeness some technical details of the LO evolution 
of the non-singlet DA are provided in Appendix lAl 



II. PARTON DISTRIBUTION AMPLITUDE 
AND LARGE iV c 



Partonic distribution amplitudes (DAs) are basic prop- 
erties of bound states of QCD. They are defined as matrix 
elements of quark bilinears between the vacuum and the 
hadronic state in question. For instance, the twist-2 DA 
of 7r a , tp( 2 \x), is given by 



(1) 



where u = 2x — 1, z is a coordinate along the light cone, 
[z,—z] denotes the gauge link operator, and /^(q 2 ) is 
the pion decay form factor, with the pion decay constant 



/ ff (0) = / = 86 MeV in the chiral limit. The DAs have 

the support x G [0, 1], normalization J„ dxip(x) = 1, and 
satisfy the crossing relation ip(x) = (p(l — x). Obviously 
the definition Q requires an identification of quark de- 
grees of freedom and also specification of a renormaliza- 
tion scale and scheme. 

On the other hand, the pion DA is related to the 
pion electromagnetic transition form factor in the pro- 
cess 7*7* — ► 7T°, or more generally, to processes with one 
pion and two (virtual or real) photons on external legs. 
This is a physical matrix element which does not depend 
on the renormalization scale and which is more suitable 
for our purposes. With the outgoing momenta and po- 
larizations of the photons denoted as q±,ei and q-z , one 
finds the amplitude 

r^ 7 » r (5i,<fe) = e^ Q/3 <eX<?2^7-7*(Q 2 ^)> ( 2 ) 

where the pion transition form factor F^ 1 * 1 * depends on 
the total virtuality, Q 2 , and the photon asymmetry, A, 



2 2 

Qi -92 



i< a < 1. 



(3) 



Equivalently, q\ = -^4^Q 2 , q\ = -^^-Q 2 . For large 
virtualities one finds the standard twist decomposition of 
the pion transition form factor 0, 



F T o 7 » 7 . (Q 2 , A) 



with 



JW(A) 



N c V X l-u 2 A 2 > 



8/A 

Nr. 





2 * ^\x)[l + u 2 A 2 ] 
ax- 



il 



[1 - u 2 A 2 }' 



(4) 

(5) 
(6) 



involving the subsequent DAs. The above re- 
sults hold modulo the logarithmic corrections incorpo- 
rated by means of the QCD evolution of the DAs, 
ip^\x) — > <f$r(x,Q) H|. At infinitely large momen- 
tum one reaches an ultraviolet fixed point behavior 
ip^ (x, oo) = 6x(l — x) , regardless of the value of the DA 
at some finite scale. It is important to realize that, at 
least in the perturbative regime, QCD evolution requires 
an identification of the powers in a twist expansion and 
the corresponding coefficients inherit a logarithmic mo- 
mentum dependence, J^ n \A) — > j( n \A, Q 2 ) provided 
their value is known at some reference scale Qq. Let us 
recall that power corrections are a high-energy manifes- 
tation of low-energy non-perturbative phenomena. 

In the large-A^ limit the vacuum sector of QCD be- 
comes a theory of infinitely many non-interacting mesons 
and glueballs (2^, 0], hence hadronic amplitudes may 
be calculated as tree-level processes, where the propaga- 
tors are saturated by infinitely many sharp meson (glue- 
ball) states. In our case the relevant diagram is shown in 
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Fig-El the pion first couples to a pair of vector mesons V p 
and V^, which then transform into photons. From sym- 
metry constraints, one vector meson has even G-parity 
(p-type) and the other one odd G-parity (w-type). Thus, 
for a massless pion we have 



F 7r o 7 , 7 »(Q 2 ,A) = 



J2 F vM) F vM)G,v p vM'<i) 

v p ,v„ 
(91* 



92), 



(7) 



where Fy p and Fy u are the current-vector meson cou- 
plings and G^VpV^ is the coupling of two vector mesons to 
the pion. At the soft photon point corresponding to the 
neutral pion decay 7r° — * 27 the chiral anomaly matching 
condition imposes the normalization 



F w o 7 , r (0,0) 



2Fy p (0)F^(0)G w y p yJO,0) 



v„v u 



4tt 2 /' 



M v B M l 



(8) 



This consistency constraint, realized in nature, can be 
always satisfied in models by an appropriate choice of 
the couplings. 

The large-Q expansion of Eq. Q would naively yield 
pure power corrections, assuming a finite number of res- 
onances ^1 ^jjl 01 ■ The situation with infinitely many 
resonances requires some care, since the coefficients of 
the twist expansion involve positive powers of the me- 
son masses and require regularization. We will illustrate 
below the situation with an explicit model. 

Another, more subtle, question regards the role of ra- 
diative corrections in the large- N c limit. This problem 
also arises in chiral quark models which are specific model 
realizations of the large- N c limit with explicit quark de- 
grees of freedom [23, HE H3, [H, ISE [M 
We take over the viewpoint adopted in our previous 
work [3^, 0, |35| , namely, we consider a situation where 
all perturbative radiative corrections are switched off. 
This corresponds to an identification of the power correc- 
tions at a given reference scale Qq. Thus, our calculations 
in the large-7V c limit determine the initial condition for 
the QCD evolution, of the form (p^ 1 (x, Qq). An obvious 
advantage of such an approach is that after the evolution 
the DAs comply to the known asymptotic QCD behavior. 
We discuss this issue in detail in Sect. IIVI 



V p ^ 



v m 



FIG. 1: The pion electromagnetic transition amplitude in the 
large-iV c limit. 



III. LARGE- iV c REGGE MODELS 

Now we proceed to the basic analysis of this paper, 
namely a study of Eq. J7J at large Q 2 . The idea is as 
follows: we adopt a model for the spectra and couplings, 
calculate the amplitude, and compare to Eq. (|5I6|I . For 
the spectra we take the radial Regge model 



My (n) = M Vm (n) = M 2 + an, 



(9) 



which assumes for simplicity the same spectra in the p 
and w channels, is well fulfilled in the experimen- 
tally explored region. We also assume that the M and 
a parameters are isospin-independent, i.e. are the same 
for the p- and w-type mesons. We investigate possible 
departures from this assumption later on. We also take 
constant, i.e. n-independent values 



Fv P 



F p = const. 



F„ 



const. 



(10) 



as follows from matching of the predictions of the ra- 
dial Regge model for the vector correlator to QCD [43| . 
Standard vector-meson dominance yields the relation 



F = F n 



N C F W . The matching yields 



a = 2ira = 



24tt 2 



-F 



(11) 



where a is the long-distance string tension. If we 
use F = 154 MeV from the p — > 2n decay [5{| we 
get -y/er = 546 MeV, while the lattice calculation of 
Ref. j5|| gives a similar order-of-magnitude estimate, 
%fo = 420 MeV. Condition (fill) originates solely from 
the asymptotic spectrum and is insensitive to the low- 
lying states, whose parameters (mass, couplings) may 
depart form the asymptotic values. 

Using the standard Feynman trick we rewrite Eq. Q 
in the form 



00 -i 

^>7*7* (Q, A ) = ^2 / dx 

n,n'=0 J ° 



(12) 



(F 2 /N c )G nn >(Q 2 ,A) 



{M 2 + anx + an'(l - x) + ±Q 2 [1 + uA}} 

where G nn > is the coupling of the pion to the n and n' 
states belonging to the p and u Regge towers. This cou- 
pling may involve diagonal (n — n') and non-diagonal 
(ri 7^ n') terms, hence we introduce n' = n + d. The 
double sum may then be transformed into 



d=0 m=0 



^ ] A nn ' — ^ ] ^ ] ( 1 j (A m ,m+d + ^m+<i,m)) 

n,n'=0 



(13) 



where the factor (1 — <5od/2) is introduced for the correct 
counting the the diagonal and non-diagonal terms. On 
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purely physical grounds we expect that the non-diagonal 
couplings are suppressed, because in that case the radial 
wave functions involve different number of nodes and the 
overlap is reduced. Thus it is reasonable to assume that 
G n ,n+d is decreasing sufficiently fast with increasing d 
such that Eq. Ijl2|l makes sense. Note that reasonable 
as it is, this assumption has sound consequences, as it 
makes the sum over n and n' in Eq. l|12|l finite. With no 
suppression the sum diverges logarithmically, however, 
any power-law suppression of G n ^ n+ d with d makes the 
expression well-behaved. At asymptotic Q 2 and fixed d 
the coupling G n;n +d might a priori depend on n. We 
show, however, that this is impossible, as it would lead 
to violation of the twist expansion 10} . Indeed, assume 
that G„, n+d = G n+d . n ~ n a ^> A ^g(d,Q 2 ,A). Then, 
using the Euler-McLaurin summation formula we can we 
transform the sum into an integral and get the asymp- 
totic behavior 



F„ arr (Q,A)~Y,(l 



■><),! 



dx / dn (14) 



n a ( d -Q 2 ^g(d,Q 2 ,A) 



{M 2 + an + adx + \Q 2 [1 + uA}} 

n a( - d ^ 2 ^g(d,Q 2 ,A) 

{M 2 + an + ad(l - x) + \Q 2 \1 + uA]}" 



dx g 



[1 - uA]"- 1 + [1 + uA]' 



whereas Eqs. Il 1151 enforce the leading behavior in the 
form 1/[Q 2 (1 — u 2 A 2 )]. This implies that the only natural 
way to match to QCD (in the absence of radiative correc- 
tions) is to assume that at large Q 2 and n we have a = 
and g depending only on d. In other words, the coupling 
Gn.n+d at large Q 2 and n is independent of n, Q, and A. 
Thus for simplicity we take G n ,n+d = G n +d,n = g{d) for 
all, even low, n. The obtained behavior is reminiscent of 
the asymptotic independence of the coupling Fy of Q 2 
and n in the case of the two-point vector correlator. It is 
worth stressing that each term in Eq. (|12|) goes as 1/Q 4 
and it is the summation over n which changes the power 
to 1/Q 2 , yielding the leading twist DA properly. Any 
finite truncation does not do the job. In order to gener- 
ate the 1/Q 2 behavior one needs the proper asymptotic 
density of states, such as in the Regge model. 
With the formula for the polygamma function, 



E 



^— ' lu + vn) 2 

71=0 V ' 



V> (-) > 



(15) 



we obtain from Eq. (|12fl the equation 



F 2 
N r a- 



^(i) ( + dx + -^Qi[i + U A] 
\ a 2a 



(16) 



M 2 , 1 9r 

+ d(l - x) + —Q 2 [l 



a 

-A). 



+(A, 

At large Q 2 we use the expansion 
1 



lA] 



^\AQ 2 
which yields 



B = 



AQ 2 



\-B 



A 2 Qi 



8F 

F»o 7 . 7 .(Q,A) = — 



l- 6 -f)g{d)j\x 



(17) 



x (18) 



Q 2 (l ~u 2 A 2 ) 



[2M 2 + a{d- !)]{! + u 2 A 2 ) 
Q 4 (l - u 2 A 2 ) 2 



Comparison to Eq. (|5I6|) gives immediately the identifi- 
cation 



.(2) 



tpr>(x) = l, tp (A) (x) 
with the normalization conditions 



1. 



A 2 



2 

Nr 



24tt 2 / 



E 



N c 
Sod 



(19) 



(20) 



g{d)[2M 2 + a{d- 1)]). 



-1/2 



as required for 3-meson cou- 



Note that g{d) 

plings, and A 2 ~ N®. The sign of A 2 is formally not 
constrained. We stress that the above results hold at the 
scale Qo at which the radiative gluon corrections are not 
present. 

Let us consider in a greater detail the simplest case 
where no non-diagonal couplings are present, g(d) = 
G5 d0 . Then Eqs. ^ yield 



G 



24^ 2 / 



A- 



M 



(21) 



With M — m p = 770 MeV the parameter A 2 is positive 
when y/a > 434 MeV. At y/a = 500 MeV we obtain 
A = 439 MeV, which agrees in the order of magnitude 
with quark-model estimates 0, • 

At Q 2 — the anomaly condition JSJ enforces the re- 
lation 




<tN c 
Airf 2 



(22) 



which is very well satisfied for a — (500 MeV) 2 , when 
the above equation becomes 0.377 ~ 0.377. If we use 
the relation m 2 = 24ir 2 f 2 /N c , then Eq. becomes 
i)^ x \z) — 3/z, with z = m 2 /a ~ 0.385, which gives 
numerically the relation yfa ~ 0.64m p = 495 MeV. This 
shows that extending the assumption of constancy of the 
the meson coupling G all the way from the asymptotic 
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region down to Q 2 — preserves the constraints of the 
model. Moreover, it leads to very reasonable results and 
allows to remarkably well determine the string tension 
from the meson phenomenology using the chiral anomaly 
matching condition. 

Let us now discuss the effects of modifying the Regge 
model. Suppose the masses M and the slope parameters 
a were different in the p and lu channels. If a p then 
the expansion ll- r >l(>[) is violated. Note, however, that in 
that case we would have different asymptotic density of p 
and u> states, which is not possible in a theory with strict 
SU(2)p symmetry. On the other hand, it is possible 
to split m p and m^. In that case still ip^ — 1, but 
higher- twist DAs become ir-dependent. We know from 
experimental data that the tu-p mass splitting is tiny, so 
that dependence should not be strong. At any case, the 
result of constant in the large- N c limit (at the scale 
Qo) seems very robust. 

Several comments referring to the interpretation of our 
results are in order. We recall that the constancy of pion 
DAs has been originally obtained by the present authors 
in the Nambu-Jona-Lasinio model |33| as well as in the 
Spectral Quark Model j^H , which are particular realiza- 
tions of the large- N c limit. A key ingredient of these cal- 
culation was the correct implementation of chiral symme- 
try through the Ward-Takahashi identities. On the other 
hand, calculations based on nonlocal quark models orig- 
inally obtained bumped distributions |2£| , close in shape 
to the asymptotic forms. Later calculations with more 
careful implementation of PCAC resulted in much flat- 
ter results |3r| . with the DA remaining non-zero at the 
end-points x = 0, 1. The trend to a flat distribution can 
also be seen in transverse lattice calculations at low trans- 
verse scales 0, • The present calculation is founded on 
more general background, using only the facts that at 
large- N c we deal with a purely mesonic theory and that 
the confined meson spectrum may be described by the 
radial Regge model. 

Another comment refers to the absence of explicit 
quarks in the present approach. Interestingly, the cal- 
culation, although referring to the partonic structure, as 
seen in Eq. Q , never explicitly introduces partons of spin 
1/2. The x variable enters from the Feynman represen- 
tation of the product of two vector-meson propagators 
and is later identified with the Bjorken x by matching to 
the QCD expressions (I'lKil) . The identification is unique, 
since DAs are universal for all Q 2 and A and are only 
functions of x (and the model parameters). 



IV. QCD EVOLUTION FOR LARGE AND 
FINITE iV c 

An important point, not only for our calculation, but 
for all model calculations, is the question of the en- 
ergy scale where the obtained predictions for the DAs 
hold. Distribution amplitudes depend on the scale, 
while our result corresponds to a fixed reference scale 



Qq. The QCD evolution is crucial, eventually evolving 

(2) 

the DAs into their asymptotic forms <^as = 6a; (1 — x), 
(p 1 ^) — 30x 2 (l — x) 2 , .... At leading order the evo- 
lution for the leading-twist component is very simple. 
One method is to use the Gegenbauer moments (see 
Appendix which evolve with the evolution ratio 

[a(Q 2 )/a(Q 2 )}^ 0) / 2 ^. Since /3 = UN c /3 - 2JV//3, and 
7n° ~ (N 2 — 1)/ (2N C ) the exponent reaches a finite value 
in the large- N c limit. As expected, the quark contribu- 
tion (depending on the number of active flavors Np) is 
only 1/N C suppressed. For instance, the second Gegen- 
bauer moment behaves as 



02(g) 

a.2(Qo) 



a(Q) 
a{Qo) 



75(JV^ -1) 
12N C (11N C -2N F ) 



(23) 



From the present calculation one gets a2{Qo) = 7/18. 
For Nf = 3 the exponent changes from 50/81 at N c = 3 
to 25/44 at N c = oo, an effect at the 10% level only. 
Thus, the (perturbative) evolution is never switched off. 
Obviously, our result for the DA cannot hold at all scales. 
As we said, it refers to a particular reference scale Qq. 
It is precisely at that situation where the DAs are con- 
stant functions of x. In order to pass to other scales, the 
evolution is necessary and its effect is strong. In partic- 
ular, the evolution changes the end-point behavior near 
x = 0, 1, for instance the twist-2 component ip( 2 ' ~ x 
near x = and p^ ~ 1 — x near x = 1 (see Refs. 
for details). This effect is very important, in particular 
for the analysis of processes with one real photon, where 
A = 1 . Then the numerators in integrands of Eq. ijoKil) 
involve powers of x(l — x) and the transition form factor 
is well-behaved only when the end-point behavior of the 
DA cancels the singularity. In short, the QCD evolution 
is mandatory if we want to compare the model predictions 
to the data. 

Having said that, we note that we can compute exactly 
the unevolved large- N c pion transition form factor, which 
amounts to carrying the x integration in Eq. (|16fl . For 
the diagonal model, g(d) = Gdda, 



Fir 



■{Q,A) 



2/ 



N C AQ 2 



M 2 



(0) 



we find 
M 2 Q 2 (l 



A) 



2a 



2a 



(24) 



This result is plotted for different values of the asym- 
metry A in Fig. [2 where the need for evolution at large 
Q 2 is vividly seen. The top curve (thick line) is for the 
un-evolved large- result of Eq. (|2*4*jl with A = 1. At 
large- Q 2 we have 



Q 2 F^o 



7 7 



(Q, \A\ = 1) 



M 2 



l og ^! +V ,(o) 

a \ a 



2/ 

+(D(l/Q 2 ). (25) 

We note the logQ 2 term, whose presence can be traced 
back to the singular end-point behavior in the twist ex- 
pansion 1 ll(il) . As a result, the pure twist expansion is 
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.25 




Q 2 [GeV 2 ] 



Aqcd = 0.224 GeV) - a rather low value which sug- 
gests the usage of higher-order evolution Q or even non- 
perturbative evolution 0, @ . Our constant DA evolved 
at LO to the CLEO scale Q 2 = 5.8 GeV 2 yields the value 
Q 2 F 7 « i7r7 (Q)/(2/ w ) = 1.25 ± 0.10 H3, higher but com- 
patible within two standard deviations to the experimen- 
tal CLEO value of 0.83 ± 0.12. This actually speaks in 
favor of small NLO corrections, however work on higher- 
order evolution should definitely be pursued in the future. 
In addition, the higher twist corrections should also be 
included in such an analysis, which is nontrivial. 
Expansion of Eq. (|2"4l at low-Q 2 yields 



FIG. 2: The pion transition form factor in the large- N c Regge 
model. Solid lines from top to bottom correspond to Eq. 12 H 
at \A\ = 1, 0.95, 0.75, 0.5, and 0, respectively. The dashed 
lines indicate the Brodsky-Lepage limit of 2/ for \A\ = 1 
(upper curve) and 4//3 for \A\ = (lower curve). The CLEO 
data points are for A = 1. 



violated. The QCD evolution cures this problem, leading 
to the asymptotic behavior in accordance to the Brodsky- 
Lepage limit, indicated with the upper dashed line in 
Fig. Therefore the QCD evolution is necessary to com- 
ply to the formal limit, as well as to the experimental 
data. More pictorially, the evolution takes the tail of 
the model calculation and brings it down to the upper 
dashed line, compatible with the data. We stress again 
that the data should not be directly compared to the un- 
evolved model results plotted in the figure, yet we notice 
that at lower Q 2 , where the effects of evolution should 
be smaller, the comparison is quite reasonable. At lower 
values of \A\ the effects of evolution are not as strong as 
at \A\ — 1, moreover, the pure twist expansion (i.e. the 
expansion in powers of 1/Q 2 without logs) holds. In the 
symmetric limit of A = we have 



F^ r (Q,A = 0) 



N c a 



2M 2 + Q' 2 
2a 



(26) 



Finally, let us mention that the determination of the 
reference scale Qo can only be done at present using 
perturbative evolution. In our previous work |33| we 
used the second Gegenbauer moment d2 = 0.12 ± 0.03 
at Q — 2.4 GeV extracted from the experiment 0] 
with the assumption at- = for k > 2. This al- 
lowed us to make the LO estimate for the evolution ratio 
a{Q)/a(Q ) = 0.15 ± 0.06. In Fig. [3 we show the cor- 
responding LO Gegenbauer evolution for N c = 3 and 
N c — oo to the scale Q 2 = 5.8 GeV 2 . As we see the 
N c -effect is tiny and is in fact comparable with the un- 
certainties induced by the evolution ratio |3^]. Let us 
note that despite the perturbative nature of our evo- 
lution the similarity of these evolved results to non- 
perturbative transverse lattice calculations with a trans- 
verse lattice size of about a± ~ 0.5 — 0.7 fm is in- 
deed striking |6l||. The LO reference scale of our es- 
timate turns out to be Q = 0.320 ± 0.045 GeV (for 



FirO- r *~ f *(Q,A) 



M 2 
a 



fQ 2 



(2) 



3a T \ a J 3a 2 
f(A 2 + 3)Q^ A3) /M 2 
36a 3 V a 



M 2 
a 



(27) 



Note that the Q 2 term is independent of the asymmetry 
A, as is also apparent from Fig. |2] The corresponding- 
slope reads 



b-K = - 



1 



7r u 77* (Q, 0) dQ 



1 ^ 

2a ^)(M±) 



(28) 



Numerically, taking M = m p = 770 MeV we get the 
value K = 1-39 GeV" 2 for a = (400 MeV) 2 , K = 
1.51 GeV" 2 for a = (500 MeV) 2 , and K = 1-58 GeV -2 
for a = (600 MeV) 2 . These Regge model estimates are in 
very reasonable agreement with the experimental values 
quoted in the PDG 58]: K = (1.79 ± 0.14 ± 14)GeV" 2 
originally reported by the CELLO collaboration [l3| . 
obtained from an extrapolation from high-Q 2 data to 
low Q 2 by means of generalized vector meson domi- 
nance, W = (f -4 ± 1.3 ± 2.6)GeV" 2 given in and 
6^ = (1.4 ± 0.8 ± 1.4)GeV" 2 given in jgjj- 

Finally, let us comment on the recent findings within 
the holographic approach [52l l53| where the light-cone 
wave functions have been computed appealing to the 
AdS/CFT correspondence and the meromorpization ap- 
proach of Ref . [54j . The holographic wave- functions cor- 
respond to a spectrum which behaves linearly in the 
mass, M(n) ~ n. Based on a conformal-based map- 
ping Brodsky and Teramond [5j| get the asymptotic 

(2) 

DA ip\'(x) — &x(l — x) while the corresponding Par- 
ton Distribution Function (PDF) is f v (x) = 6x(l — x). 
If, instead, a twist-based mapping is considered |52( 

these authors get ip^\x) — (&/tt)\Jx{\ — x) and a PDF 
f v (x) = 1. According to the meromorphization approach 

of Radyushkin [54[ one gets ip^(x) — (8/ir)y/x(l — x) 

and fn(x) — 1 for scalar quarks and ipn(x) — 6x(l — x) 
, the asymptotic DA, and f v (x) = 6a; (1 — x) for spin 
1/2 quarks. This latter situation is exactly the kind of 
situation that was found in the NJL quantized on the 
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1.5 



■■ 5.8 GeV, 



=3 , 

N c =oo , = 5.8 GeV 
6x(1 -x) (Asymptotic Result) 
(f % (x, Q ) = 1 (Initial condition) 




FIG. 3: The pion distribution amplitude evolved to the scale 
Q 2 = (2.4GeV) 2 for the cases N c = 3 and iV c = oo. The 
value for the evolution ratio a(Q)/a(Qo) = 0.15 is based on 
the analysis of the CLEO data of Ref. [Tsj . We also show 
the un-evolved DA, <p n (x,Qo) = 1, and the asymptotic DA, 
ip*(x, oo) = 6x(l — x). 



light cone 27]. An asymptotic PDA suggests a refer- 
ence scale Qo = oo or the lack of LO evolution of <fi 2 , 
which assumes the asymptotic form at all scales. On the 
other hand the corresponding PDF yields a 100% mo- 
mentum fraction carried by the quarks, while for Qo = oo 
one expects that all momentum fraction is carried by the 
gluons. This poses an interpretation difficulty for these 
approaches which should be cleared out. 



V. CONCLUSION 

We summarize our main points. The radial Regge 
model predictions for the 7r°7*7* amplitude can be 
matched to the QCD twist expansion in the absence of 
radiative corrections in the large- N c limit. The match- 
ing requires that the coupling of the pion to a pair of 
Regge p and u> mesons, G^vv, is constant for highly 
excited Regge states and large momenta. The twist-2 
pion distribution amplitude is then found to be constant 
in the cc-variable, conforming to the earlier predictions 
made in a class of chiral quark models. Thus, we pro- 
vide an explicit example of quark-hadron duality for an 
exclusive process. Moreover, higher- twist DAs may or 
may not be constant, depending on the details of the 
Regge model. If m p — m^, then the twist-4 DA is also 
constant. The Regge model couplings are constrained 
by matching the chiral anomaly for real photons and the 
high-momentum behavior, cf. Eqs. (|2t)l21(l for highly vir- 
tual photons. If we further assume a model with constant 
diagonal couplings, a consistency relation (|22f> involving 
the string tension is found, which is well supported by 
the data. As a result the string tension is found to scale 
with the square of the p-meson mass. The Regge model 
results is very reasonable predictions for the slope of the 
pion electromagnetic transition form factor at zero mo- 
mentum. The estimates depend weakly on the string 



tension, and compare quite well to the existent measure- 
ments. Finally, we have noted that similarly to the case 
of chiral quark models, the QCD evolution is necessary 
in the present Regge model to achieve the correct large- 
momentum behavior of the pion transition form factor. 
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APPENDIX A: LO EVOLUTION OF DA 

The LO-evolved distribution amplitude reads P, 0] 

oo ' 

<pW{x,Q) = 6x(l-x)J2c 3 J 2 (2x-l)a n (Q)lAl) 

where the prime indicates summation over even values of 
ri only. The matrix elements, a n (Q), are the Gegenbauer 
moments given by 



a n (Q) 



3/2 

with C n denoting the Gegenbauer polynomials, and the 
non-singlet anomalous dimension reads 



2n + 3 



a(Q) 

3(n+l)(n + 2) \a{Q ) 



7i°V(2/3o) 



dxC^ 2 (2x-l)ip^(x,Q ), 



(A2) 



ri 0) = -2C f 



n+l 



(n + l)(n + 2) 



_ i 



k 

fe=i J 



(A3) 



with C A = N c , C F = (N 2 - 1)/(2N C ), and N F being the 
number of active flavors, which we take equal to three. 
The LO running is then a(Q) = (4tt//3 )/ log(Q 2 /A^ CD ). 
Taking as initial condition 



^ 2 Hx,Q ) = l 
one we gets immediately 



a n {Q) 



2ri + 3 



am 

3(n+l)(n + 2) \a(Q ) 



7™/(2/3o) 



(A4) 



.(A5) 
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